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Tomographic techniques have great abilities to interrogate the combustion processes, especially when it is com- 


bined with the physical models of the combustion itself. In this study, a data fusion algorithm is developed to in- 


vestigate the flame distribution of a swirl-induced environmental (EV) burner, a new type of burner for low NOx 


combustion. An electric capacitance tomography (ECT) system is used to acquire 3D flame images and computa- 


tional fluid dynamics (CFD) is applied to calculate an initial distribution of the temperature profile for the EV 


burner. Experiments were also carried out to visualize flames at a series of locations above the burner. While the 


ECT images essentially agree with the CFD temperature distribution, discrepancies exist at a certain height. When 


data fusion is applied, the discrepancy is visibly reduced and the ECT images are improved. The methods used in 


this study can lead to a new route where combustion visualization can be much improved and applied to clean 


energy conversion and new burner development. 
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Introduction 


As the requirements for pollution control and energy 
saving have become more stringent in recent years, the 
study of low-emission and high-efficiency burners is in- 
creasingly necessary. It is believed that low-heat value 
gaseous fuels, such as blast furnace gas for gas steam 
combined cycle power generation, not only can solve the 
problem of electricity production in the iron and steel 
industry, but also make breakthroughs in energy saving, 
emission reduction, and environmental improvement [1]. 

New burners for low-heat value gaseous fuels have 
been studied for many years. In the 1990s, ABB solved 
the problem of combustion stability with a complex con- 
trol system for a low-heat value gaseous fuel burner [2]. 

Among many new technologies, the swirl-induced en- 
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vironmental (EV) burner has shown favorable features 
such as low emission generation and good combustion 
stability. However, as a relatively new type of burner, 
detailed investigations are still needed to understand the 
mechanisms of the combustion process. The spatial dis- 
tribution of the flames in the EV burner is of great im- 
portance in assessing the combustion conditions, as it is 
for common combustors. As combustion very often oc- 
curs with variations or fluctuations in parameters such as 
combustion intensity in various locations, monitoring of 
flames should not only be of spatial distributions, but 
also be able to follow the variations with a certain tem- 
poral resolution. 

To achieve this goal, various techniques are sought to 
visualize flame distributions in combustion chambers. 
This will be the first and necessary step leading to more 
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difficult tasks such as detecting the locations of poor 
combustion zones that tend to reduce the combustion 
efficiency. It would also highlight the locations with ex- 
cessively high temperature that would increase nitric 
oxide emissions (NOx) [3]. Moreover, probing details of 
combustion should require observation of the internal 
structures of the flames. As a relatively new technology, 
electric capacitance tomography (ECT) is an attractive 
approach for its ability to map the entire cross section of 


a combustion zone with a very fast data acquisition speed. 


ECT is by nature sensitive to electrical properties of the 
flame, and maps permittivity of the flame that is closely 
related to the intensity of the combustion. 

Electric capacitance tomography, ECT for short, is a 
recently developed non-intrusive visualization technolo- 
gy. With distinct advantages, such as non-intrusive sens- 
ing and fast data acquisition speed, ECT reconstructs the 
permittivity distribution in the measuring zone using cer- 
tain algorithms on the measured capacitance data. The 
ECT method has been intensively studied and is one of 
the most promising visualization methods for parametric 
measurements in multiphase systems. Especially, Water- 
fall et al. [4], He et al. [5], Liu et al. [6] apply the ECT 
method into the visualization of the spatial distribution of 
a flame, which provides an effective approach for under- 
standing the underlying physical and chemical mechan- 
isms of the combustion process. 

Depending on the materials being measured and the 
working environment, ECT can use various types of 
sensors. Figure 1 illustrates an ECT system with a cylin- 
drical sensor. A data acquisition device collects data from 
the sensor and sends them to a PC, where the data are 
analyzed and images are reconstructed. 


electrode reconstruction image 


Insulated aS 


layer 


medium 
Imaging system 


Sensor system 


Fig. 1 A sketch of typical ECT system 


According to current theories, flame ionization con- 
sists of thermal ionization and chemical ionization. 
Chemical reactions and high temperature may excite 
substances in the flame and generate electrons, ions and 
electric dipoles. Thermal ionization is determined by 
flame temperature and ionization potential. On the other 
hand, chemical ionization is the process that is closely 
related to endothermic and exothermic reactions of the 
combustion. For a methane flame, for example, the posi- 
tive ions in high concentration are C3H;° , H;O*, CHO’," 
H;O,", NO’," C,0>H’, HO’, H;0,", etc. , however, the 
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negative particles in the flame are mainly electrons and a 
few negative ions. 

As far as 3D ECT techniques are concerned, there 
have been active investigations in both sensor design and 
image reconstruction algorithms recently. While a major- 
ity of the efforts have focused on 3D multi-phase imag- 
ing, visualization of the details of a combustion process 
and flame characteristics still needs substantial develop- 
ment. Ye [7] investigated volumetric and planar ECT. In 
his study, a planar sensor was built with a matrix of three 
by four rectangular electrodes laid on a plastic base. The 
Landweber iterative algorithm was used for the image 
reconstruction. Images of test objects, such as wooden 
blocks, were reasonably well reconstructed. Depth of 
detection was also observed and found to be at least 8 cm 
for a sensor size of 17 by 17 cm. Soleimani et. al. intro- 
duced an ECT sensor of 32-eletrodes arranged in four 
planes with eight electrodes on each plane [8]. They im- 
plemented a temporal correlation between 3D ECT im- 
ages to form 4D maps. 3D ECT has also been used for 
real time multiphase flow visualization. Warsito and Fan 
[9] applied the 3D ECT technique to gas-liquid-solid 
fluidized bed measurements, and successfully revealed 
the flow of bubbles in a multiphase system. As for flame 
imaging, Gut and Wolanski carried out a preliminary 
study on 3D flame imaging using a cylindrical ECT sen- 
sor of 18 electrodes distributed in three planes [10]. This 
article showed a group of flame images in 3D shapes. 
However, the inner structure of the flames was not re- 
ported in the measurement results. All these works have 
shown certain fundamental features of 3D ECT. 

Previous work of our study [11] has tested a planar 
circular sensor that serves well the purpose of 3D flame 
imaging, being capable of revealing different zones of a 
flame. Being a preliminary study, the full characteristics 
of such a sensor have not yet been fully examined. Fur- 
ther work is thus expected. 

The images reconstructed by ECT are presumed to be 
related to the combustion intensity that is in one way or 
another related to the temperature distribution. Based on 
the above assumptions and to improve the ECT technique 
to reflect the temperature profile more rigorously, it will 
be beneficial to integrate the information of the combus- 
tion itself. For this purpose, in this investigation a nu- 
merical simulation of combustion is performed and its 
data are used as a source for data fusion with the ECT 
measurements. Details will be described in the following 
sections. 


Models and methods 


CFD model 


The objective of this study is to monitor combustion 
or the flames in an EV burner by means of a non-intru- 
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sive visualization technique, namely ECT. Basically an 
EV burner is a conical shaped burner, with two air slits 
opened at opposite sides of the cone wall and the com- 
bustion air flows tangentially from the two slits into the 
cone of the burner. In the middle part of the slits there are 
small holes acting as the inlet of the gaseous fuel. The 
fuel and combustion air should mix well for the complete 
combustion of the fuel. Due to the varying radius of the 
cone, the swirling gas will form a stable recirculation 
zone and a stable combustion flame in front of the recir- 
culation zone, see depiction in Fig.2. 


(a) External air inlet 


(b) Combustion chamber 


Fig. 2 The structure of EV burner 


A Computational fluid dynamics (CFD) simulation is 
carried out to study the combustion process in the EV 
burner. In the current work, the Reynolds averaged Navi- 
er-Stokes equations combined with the energy equation 
are solved to obtain the flow field and the temperature 
profiles in the combustor. The continuity, momentum and 
energy equations are listed below: 


V-(p7)=0 (1) 


V-(pi¥) =-Vp +V-(r) + p8 (2) 


ve =p) =v ievr -Zn (E) A 
i 

where p is density, v is the time average velocity, T 
denotes the stress tensor and g is the gravitational ac- 
celeration. For the term E = h — p/p+v’/2 in equation 3, h 
represents the sensible enthalpy and p is the mean pres- 
sure. kepis the effective conductivity and J represents 


the diffusion flux of species. 

Very often the standard k-e turbulence model is ap- 
plied to close the Navier-Stokes equations. The transport 
equations for k and ¢ are given below: 


v (ote) =v (urt )ve]+a +O, wy 


E e 
+ Cy, re, + C3.G5)— Cre T (5) 


where k and ¢ are the turbulent kinetic energy and the 
turbulent dissipation rate respectively. u is the air viscos- 
ity and u, represents the turbulent viscosity. o, = 1.0 and 
o= 1.3 represent the turbulent Prandtl numbers for k and 
€ respectively. G; and G, denote the generations of tur- 
bulence kinetic energy due to mean velocity gradients 
and buoyancy respectively. The constants of C;, and Co; 
are set at the default values 1.44 and 1.92. However, oth- 
er turbulent models, such as the Reynolds stress model, 
can also be considered. 

In order to calculate the turbulent flame in the com- 
bustor, the non-premixed model is adopted and the go- 
verning equations are shown in equations 6 and 7. In the 


equations f is the Favre mean mixture fraction and 


f? is mixture fraction variance. The constants of op C, 
and C4are defined as 0.85, 2.86 and 2.0 respectively. 


V-(pv f)=V ‘aaa 
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(6) 
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The commercial CFD software Fluent is used to solve 
the equations. The SIMPLE algorithm is chosen to 
couple the pressure and the velocity fields and 
second-order schemes are used to solve the discretized 
form of the above equations. To deal with the high tem- 
perature inside the combustor, the Discrete Ordinates 
model is adopted to consider the radiation effect. The 
inlets of the fuel (CH,) and air are treated as the velocity 
boundary and the outlet of the computational domain is a 


pressure outlet. In addition, the variable f is set to 1.0 at 


the fuel inlet and 0 at the air inlet, respectively. Also, the 
same combustion conditions, e.g. structure of the burner, 
air feed rate and fuel feed rate, are applied to both the 
CFD simulation and the ECT measurements to acquire 
comparable data. 

The temperature distribution calculated by CFD could 
be regarded as a good reference for the ECT measure- 
ment, either as prior knowledge for the ECT image re- 
construction, or as a comparison to test and verify the 
accuracy of the reconstruction algorithm. In this way, 
when the parameter distributions from the CFD simula- 
tions are successfully integrated into the images from the 
ECT experiments, better results should be expected. 


The ECT Sensor 


A planar circular sensor with eight electrodes was de- 
signed and fabricated to detect the material distribution in 
a 3D space, and to be suitable for a burner, as shown in 
figure 3. In the central part of the base is a hole to host 
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1, electrodes; 2, stainless steel base; 3, crystal glass; 4, data line; 5, supporting flange; 6, shielding; 7, EV burner 


Fig.3 Photos of the eight-electrode ECT sensor 


the EV burner. The variation of the distribution of the 
materials above the electrodes will cause the corres- 
ponding capacitance variation that will be measured by 
each pair of electrodes. Therefore, the 3D distribution of 
materials can be detected. For the same reason, when 
there is a flame, the permittivity distribution above the 
electrodes will also be sensed by the electrodes. 

The sensor in figure 3 has eight electrodes evenly 
spaced on the top plane, beneath the electrodes is an in- 
sulation plate made of crystal glass, further down is a 
stainless steel base. The EV burner is installed through 
the central hole. 


3D Sensitivity map 


A sensitivity map is essential for ECT image recon- 
struction. In this study, a two-step numerical procedure is 
used to generate the sensitivity map: 

a. The electrodes are excited one by one. For each ex- 
cited electrode, the electric field is calculated based on 
electromagnetic models via the finite element method. 
This process continues until every electrode is excited 
and eight electric fields are obtained corresponding to the 
eight electrodes. 

b. For each pair of the above calculated electric fields, 
the inner product is performed to obtain the sensitivity 
map for each voxel cantered at (x;, Yi, Zi). 


SCi j,xX,y¥,Z) = E; x(x, y,Z) E; (x,y,z) 
+E; (x,y,z) E; (x,y,z) (8) 
+ E; (x, y,z): E(x, y,Z) 
where S(i_j, x, y, Z) stands for the element in the sensi- 
tivity map corresponding to the electrode pair i, j, and a 
voxel location (x; yi, zi); Ex (x, y, Z) is the electric field 
matrix at (x, y, Z), with the first subscript indicating the 


number of the electrode and the second subscript indi- 
cating the component of the electric field. 


Image reconstruction and fusion method 


The ECT image reconstruction procedure consists of 
two processes: the forward problem and the inverse 
problem. In a mathematical notation, when signal noise is 


ignored, the common ECT imaging model can be ap- 
proximated as: 
SG = C (9) 
where S stands for the sensitivity map, or measure- 
ment matrix; G stands for the permittivity distribution 
vector, also called the image that contains the grey values 
for each voxel; C holds the capacitance values, measured 
or simulated. 

When reconstructing the images, an iterative proce- 
dure can offer more control over the quality of the image 
and the magnitude of the error. In this paper, the Land- 
weber algorithm is employed to reconstruct the images. 
The general equation describing the Landweber algo- 
rithm is: 

G; =G; +aS!(C-SG;) (10) 


Where æ is the step length of the iteration used to con- 
trol the convergence speed. When the initial value Go, 1.e. 
the first image in the iteration, is not known it can often 
be assigned with random values or simply zeros and the 
iteration process then converges to the value of the min- 
imum norm square. This may cause bigger errors or a 
longer iteration time. 

An effective way to judge the convergence for the 
ECT inverse problem is to set a tolerance error value & 


for the iteration: 
Gii-G;< č (11) 

When the condition of equation 11 is fulfilled, the ite- 
ration process is completed. 

For the purpose of data fusion, there can be a variety 
of ways of introducing CFD data into the image recon- 
struction process. In this study we propose using CFD infor- 
mation as prior knowledge and assigning CFD data to Go 

Go(x, y, z) = CFD_data(x, y, z) (12) 
where CFD_data are in a normalized vector contain- 
ing the temperature data for each voxel. 


Experimental 


Experiments to image the flames were carried out and 
figure 4 shows the experimental apparatus. 
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1. Methane cylinder, 2. Air cylinder, 3. Pressure gauge, 4. Pressure regulator, 5. Flow rate valve, 
6. Float meters, 7. ECT sensor, 8. Combustion chamber, 9. ECT data logging device, 10. Computer. 


Fig. 4 Experimental apparatus 


Table 1 Parameters of the Cyclone 


Number 2 
o . Length/ mm 70 
Air-intake slit 
Width /nm 2 
Intake area/ mm? 280 
: Number 2 
Fuel gas inlet 
Length/ mm 70 
Number 8 
Gas hole Radius/mm 0.5 
Exit area/mm? 6.28 
Angle/* 34 
Half-angle of the Upper end diameter /mm 60 
combustion cone Lower end diameter /mm 8.04 
Height/mm 65 
Inner diameter /mm 185 
flame tube p 
Height/mm 200 


In this study an EV burner is designed and the specific 
parameters are shown in table.1 The fuel is high purity 
methane (volume fraction, 99.95%) and dry air is the 
oxidant. The feed rates of the fuel and air are controlled 
by the valves in the flowmeters. The ECT sensor is fitted 
onto the gas burner to form an integrated unit, shown in 
figure 4. For the eight-electrode ECT sensor, the outer 
diameter of the electrode array is 124mm, with an inner 
diameter of 74mm. The electrodes are made of stainless 
steel with a 4 mm thickness and mounted on a crystal 
glass plate of 12 mm thickness. Below the glass plate is a 
metal plate serving as a shield to prevent interference 
from the outside environment. Below the metal plate 
there is another crystal plate to insulate the electrode 
wiring from the metal plate. For further reduction of en- 
vironment interference, a cylindrical shielding was also 
used to isolate the combustion space from the outside 


space. A range of methane fuel feed rates were tested, 1.e. 
from 30 L/h to 180 L/h; and feed rates of the combustion 
air ranges from 500 L/h to 3000 L/h. 

To avoid signal saturation during the measurement, 
calibrations are conducted before the experiments con- 
sisting of a low range and a high range measurement. The 
low range calibration is a normal measurement with an 
empty combustion chamber. The high range calibration is 
performed using sand with a series of moisture content 
that yields a range of permittivities from 3 to 15. After a 
number of trial tests, it was found that for the combustion 
conditions available, sand with a permittivity of above 8 
for the high range calibration was sufficiently high to 
prevent signal saturation. 

In this study, 1000 frames at a data sampling rate of 40 
frames per second were normally captured for each expe- 
rimental run, which is enough to provide adequate tem- 
poral resolution. With the measured capacitance and the 
3D sensitivity map, a series of images of the flames are 
reconstructed. 


Experimental results 


A set of typical experimental results are presented here 
when the feed rate of methane is 90 L/h and the feed rate 
of the combustion air is 1000 L/h, i.e. the fuel/air ratio is 
1:11. For this case, the 3D images were reconstructed 
using the 3D sensitivity map and the Landweber algo- 
rithm. To provide clear views of the inner structure of the 
flames, a series of slices of a typical 3D image are pre- 
sented in figure 5, each slice corresponding to a different 
cross section along the axis of the combustion chamber, 
at different heights from the surface of the electrodes of 
the ECT sensor. There are three typical cases in figure 5, 
with similar characteristics. 
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H=20 mm H=60 mm H = 100 mm H= 140 mm H=180 mm 
1 
2 
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Fig.5 Slices of the 3-D images of the inner structure of a flame at different heights H in different time 


The basic features of figure 5 are clear. First, for low 
positions (e.g. H = 20 — 60 mm), the high permittivity 
(high grey values, shown as red) parts are spread out 
above the sensor. Second, at a certain height (e.g. H is 
100 mm), the image is distributed in a circular manner. 
Third, the grey values in the cross sections are not uni- 
form, with strong and weak portions, which implies an 
asymmetrical intensity of the flame. This can be consi- 
dered reasonable as the combustion is of a turbulent na- 
ture. Fourth, at high positions (e.g. H = 140 — 180 mm), 
the circular structure reduces and eventually merges into 
a small region of similar intensity, which may indicate 
the top region of the inner cone of a flame. On the whole, 
the images display the spatial distributions of the flames 
quite well. 


CFD results and the effect of data fusion 


CFD simulations were conducted by solving the con- 
trol equations for the Navier-Stokes equations. For data 
fusion purposes the simulations adopted the same condi- 
tions as those used for the ECT measurements. The in- 
ternal flow field of the burner is based on the pressure 
solver and the pressure-velocity coupling algorithm. The 
standard wall function is used near the wall surface. The 
intake velocities include the velocities of the air and the 
fuel, and the inlet temperature is set to 300K. The burner 
exit condition is set to the flow boundary. All the wall 


surfaces are assigned the solid walls. An upwind differ- 
ence scheme of second order accuracy is adopted. Mul- 
tiple grids are used to accelerate the convergence of the 
computation process. 

The temperature distributions at a series of cross sec- 
tions are presented in figure 6 below. For clearer visual 
effects, different scales are used in figure 6, i.e., from 
500K to 1550K for the top row and from 600K to 1200K 
for the bottom row. 

The distribution of the temperature in the combustion 
zone will reflect the flow of the combustion media and 
the products, the shape, the location, and many other 
important characteristics of the flame. As is shown in 
figure 6, the flow in the burner appears to be rotating. In 
a certain range of heights, i.e. H from 10 mm to 30 mm, 
there is a band of high temperature region spreading from 
the central zone towards the edge. This is characteristi- 
cally in agreement with the ECT images for those heights. 
However, the CFD data more clearly displayed a swirl 
like shape at H = 10 mm and at H = 30 — 180 mm, which 
very likely reflects the effect of the two air inlet slits at 
the burner wall from which the fuel and air enter the 
burner in a tangential direction. At some intermediate 
heights (e.g. H = 60 — 100 mm), the temperature appears 
to have a circular shape, which resembles well the ECT 
image for H = 100 mm. After a certain height the tem- 
perature tends to decrease and be more uniform, as at H = 
140 mm and 180 mm. 
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By comparing figure 5 with figure 6, one can see that 
the images by the two methods roughly agree for certain 
heights, say 20 mm — 100 mm, while the results for H = 
140 mm and 180 mm show certain discrepancy: the high 
intensity areas in the ECT images are small but the CFD 
profiles are largely uniform. The reason might be that in 
the higher positions the sensitivity of ECT is weak, thus 
the ECT images show small area of high temperature. 
While the CFD results are not related to the sensitivity of 
ECT detection, they may be more faithful to the real 
condition and display larger areas of more evenly distri- 
buted temperature profiles. Based on this consideration, 
one may expect improved ECT results by data fusion 
with the CFD data. 

Data fusion was performed using equation 12, where 
the CFD data provide the prior information as Go, and the 
iterations were performed using equation 10. Usually Gy 
was assigned random values or the values of a linear 
back projection of the measurement. When Gy is decided 
by data fusion with complementary data, i.e. CFD data in 
this study, the results should be improved. Figure 7 
shows the images after data fusion. Improvement can be 
seen for the ECT images, e.g. for H = 60 mm and 180 


H=10mm H=20 mm 


H=60 mm = 100mm 


Fig. 6 
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mm, where clearer circular distribution is displayed for 
the former and the image is more uniform for the latter. 
Therefore, by applying data fusion, the ECT images come 
to a higher degree reflect the temperature distributions. 


Conclusions and Discussions 


In this study a planar circular ECT sensor was devel- 
oped to be used with tomographic techniques for EV 
burners aimed at 3D visualization of flame structures, 
while a data fusion algorithm was derived to improve the 
results. Experiments were conducted to obtain ECT im- 
ages of the flames, while the temperature distribution of 
the flame under the same experimental conditions was 
calculated using CFD simulations. ECT images clearly 
revealed the structure of the flames. Albeit with certain 
discrepancies, in general good agreement was found be- 
tween the ECT results and the CFD data. In addition, 
data fusion effectively improved the ECT images and 
reduced the discrepancies between ECT and CFD results. 

Future study will be needed to find more optimized 
sensor structures and to develop more accurate images 
fusion algorithms for better results. 


H = 30 mm H = 40 mm 


= 140 mm = 180 mm 


Slices of the temperature distribution of the flame by CFD 


=20 mm H = 60 mm 


H = 100 mm 


= 140 mm = 180 mm 


Fig.7 Fusion images of ECT and CFD data 
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